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Abstract 

It is demonstrated that w _4//3 law of the power spectrum with the angular 

frequency u in dissipative flows is produced by the emission of dispersive waves 

from the antikink of an congested domain. The analytic theory predicts the 

spectrum is proportional to lu~ 2 for relatively low frequency and u; -4//3 for 

high frequency. 
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Recently, much attention has been attracted to collective dynamics of dissipative particles 
[1,2]. In particular, physics of granular flows [3-6] and traffic flows [7,8] are developing 
subjects. In such dissipative flows, we often observe the coexistence of congested regions 
and dilute regions. It is important to know the mechanism of the emergence of congestion 
of traffic and granular flows. Although we have some exact results on the formation of 
congested domains in one-dimensional traffic flow [9,10], we still do not understand the 
details of the fluctuation of dissipative flows. 

In experiments of dissipative flows, we usually measure the power spectrum which is 
the Fourier transform of the auto-correlation function. It is known that traffic flows and 
granular flows in a pipe have the power spectra obeying law with the angular frequency 
lu [11,12]. Several years ago, Moriyama et al. [13] have confirmed that granular flow in a 
pipe should have the spectrum with (3 = 4/3. We also expect that the power spectrum 
obeying to~ 4 ^ 3 law is universal for dissipative flows in the coexistence of congested-flow and 
sparse-flow. [3,13-16] This law is robust in the experiments of granular flows, which can be 
observed without tuning of a suitable set of parameters. [13,16] 

Although the previous papers [3,13] proposed the mechanism of tu~^ 3 law, their deriva- 
tion might be incomplete. We can list several defects in their derivation: (i) They assumed 
that the system is in a weakly stable region of homogeneous state. However, cu~ 4 ^ 3 can be 
commonly observed in the case of the coexistence between congestion and the sparse-flow. 
The power exponent (5 is drastically small when there are no definite domains in systems. 
[15,16] (ii) The experiments [13,16] suggest that co~^ 3 law is robust without fine-tuning 
when phase separations take place, but the theory assumes that the system is in the vicinity 
of the neutral curve of the linear stability analysis, (iii) The theoretical spectrum depends 
on the wave number but there is no wave-number dependence in the actual observation in 
experiments. [3,13] (iv) Although the theory assumes that the relaxation process of internal 
structures, it is not clear what the relevant relaxation process is. Therefore, one is skeptical 
of the validity of the previous theory to explain cj~ 4 / 3 law. 

Recently, Takesue et al. [17] have solved a kink-diffusion problem in the totally asym- 
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metric simple exclusion process (TASEP) [18] and derived uj~ 3 ^ 2 law of the power spectrum. 
Although TASEP contains only a kink which connects one congested domain with a di- 
lute region, their analysis is suggestive to understand more realistic situations in traffic and 
granular flows. In this Communication, we thus try to re-derive tu~ 4 ^ 3 law in the case of 
coexistence between congestion and sparse-flow. 

In order to proceed the analysis we should recall that all of one-dimensional models for 
traffic and granular flows in weakly unstable regions can be described by trains of quasi- 
solitons stabilized by small dissipations. [3,19-21] In general, a dilute region is connected 
with a congested region by asymmetric interfaces [3,20,21] which may be characterized by 
the soliton equation. [19] We call a front interface the kink and a backward interface the 
antikink. The antikink is not stable in the actual situations and emits dispersive waves 
backward. The waves are caught by the next domain. In the simplest situation, we can 
ignore the widths of the kinks and antikinks which may be much smaller than the typical 
domain size. 

From the observation of experiments the power spectrum may not be related to the 
formation process of domains but be characterized by the emission of dispersive waves from 
an antikink. Thus, we ignore the formation of a congested domain but focus on the decay 
process of the domain. We also map the model onto a one-dimensional space, where the 
position fixed in an experimental system is denoted by x and the system size is L and the 
boundaries are located at x — ±L/2. For simplicity, we place a detector to measure the 
power spectrum at x — 0, i.e. the center of the system. Let us introduce the packing fraction 
(f>(x,t) = n(x,t)/no where no is the maximum density. 

If we assume that an idealistic congested domain exists in the system at time t — 0, the 
packing fraction is given by t = 0) = 1 between x = xq and x = xq + /, and <f>(x, 0) = 
for otherwise, where / and xo are the size of the domain and the position of an antikink at 
t — 0, respectively. The equivalent expression is 
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On the other hand, the antikink is unstable because of the dispersion of propagating velocity, 
though we can ignore such the effect for the kink. Thus, we assume that the time dependence 
of 0(0, t) can be described by 
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where c is the average speed of domains and £ is a characteristic length scale of the dispersion 
relation. We note that the shape of domain at time t is no longer idealistic one but is decayed. 

With the aid of Wiener-Khinchin theorem, the power spectrum I(u) and the auto- 
correlation function C(t) can be written as 

I roo 

I{u) = -= / dtc^C(t), C{t) =< 0(0, 0)0(0, t) >, (3) 

where the ensemble average in eq.(3) is interpreted as the average by the position of the 
antikink xo- Because the domain propagates with the constant speed Co if we neglect the 
dispersion, the existence probability of domains should be uniform except for the boundary 
regions. Thus, we may assume the probability distribution function P(x ) = 1/L and 
C(t) = lf^ L 2 /2 dx (f>(0, 0)0(0, *). Substituting eqs.(l) and (2) into (3) we obtain 
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Here, J (t) in eq.(4) can be calculated as 



where we use the formula Y^=\ cos nx/n 2 = ir 2 /Q — nx/2 + x 2 /A. Thus, Iq(u) = 
^ X!^ dte iujt {l 2 / L 2 + J (t)) becomes 

W = v^±V) + ^c-. (9) 

The evaluations of J\(t) and J2W are nontrivial. When we assume c^t <C L the summa- 
tion in J±(t) can be replaced by the integral. From the expansion by c t/£ we obtain 

2/3 



sinz. 



' -(£ 2 c t) 1/3 + ^, (10) 



3Lr(4/3) vs u ' 6L' 

where we use / °°ds(l - cos^)/^ 4 / 3 = ?r/r(4/3) and $™dzwa.z/z = n/2 with the Gamma 
function r(z). The corresponding Fourier transform of Ji(t) is thus given by 

On the other hand, for / ^> c t, can be evaluated as 
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Here J21XO and J22(£) are respectively given by [22] 
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with b = I; 2 c /I 3 , and 
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where i-F 4 [ai; b u b 2 , h, 6 4 ; z] = T,kL ( a i)kZ k /[(h)*; ' ' ' (h)kk\] with (a) k = a (a + 1) • • • (a + 
fc — 1) is the generalized hypergeometric function. 

The expressions of J 21(f) and J 22(f) are complicated. For example, the Fourier transforms 
of J 2 i (t) is given by [22] 
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It should be noted that the Fourier transform of J -off) becomes complex. This is because 
the expression we obtain depends on the initial condition and the choice of the frame. If 
the system is Galilei invariant, such the term should be zero. Therefore, we regard the 
contribution from 322(f) as unimportant. 
In the limit of uj — > 0, Ii(w) becomes 

I2H * ^hfu,) -^(^co) 1 ^" 4 / 3 . (16) 

It is notable that this asymptotic expression of 12(00) is canceled with the term proportional 
to c<j _4 / 3 in h(oj). That is, the spectrum obeying uj~^ 3 disappears and I(uo) ~ uo~ 2 in the 
limit of uj — > 0. 

On the other hand, though li\(u)) is singular in the limit of uo — > 00, l2i(w) is regular for 
enough large to. In fact, one can obtain the analytic expansion of Ji\(f) near bt = 0.001 as 
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2J 2 i(t)/ir ~ l + a(bt- 0.001) +O((6t-0.001) 2 ) ~ 1 + abt + ■ ■ ■ with a = 0.000229538. If we 
replace J2i(i) by this approximate function, we obtain the approximate Fourier transform 
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for large uj. 

Thus, we obtain the power spectrum I{uj) = Io(u) + h(oo) + as 

for uj 7^ 0. For large cu, I(u) is dominated by the term proportional to u^^ 3 as 

Thus, we derive the spectrum obeying lu~^ 3 . Figure 2 shows the comparison of eq.(18) with 
eq.(19), where we can see the tail obeying cj -4 / 3 for large u, while eq.(18) seems to obey u~ 2 
for small uj. It is obvious that both expressions (18) and (19) become identical for larger uj. 

It should be noted that the actual process includes many other factors for larger u and 
smaller uj. In experiments, decays exponentially for larger uj, because the initial state 
is not in an idealistic domain as we have assumed here. To reproduce the full shape of 
spectrum we need to contain the formation process of domains for our analysis. That is an 
important future problem to be solved. Nevertheless, we believe that our picture presented 
here captures the essence of physics and clarify the mechanism of emergence of u~^ 3 law. 

In this Communication, we have demonstrated that the main process to produce uj~ a I 3 
law is the emission of the dispersive wave from an antikink. This result is universal when 
isolated congested domains exist in a dissipative flow. Through the analysis, we have revised 
the previous uncertain picture. 
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FIGURES 
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FIG. 1. A schematic picture of the propagation of a domain (its size I) with the speed cq. The 
dispersive wave emits from the antikink interface at Xq. 
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FIG. 2. Log-log plots of eqs.(18) and (19) as the frequency spectra. We adopt the parameters 
l/L = i/L = 0.3 and Lb/c Q = L/l = 10/3. . 
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